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It is shown that the observability of a large class of operations on mixed states is fundamentally 
limited. We consider trace preserving, unital operations. This class includes unitary and perfect 
premeasurement operations. An upper bound on the trace distance between an untransformed state 
and a state transformed by one of these operations is derived. The bound is dependent only on the 
purity of the state. In the case of maximal mixedness, the bound implies all operations of this class 
are unobservable. 
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I. INTRODUCTION 

Given a particular quantum state that is subjected to 
a class of operations on the state, does the state change? 
And if so, how observable is that change? This question 
is related to the study of the preservation of informa- 
tion subject to quantum processes [Tj-Q by Nielsen et al. 
Their work focuses on the slightly different question of 
quantifying changes to quantum states given specific op- 
erations. By asking our initial question instead of that 
raised by Nielsen et al. the properties of the state are em- 
phasized. We shall in particular explore the relationship 
between structural properties of a state (i.e., entangle- 
ment) and its dynamics. 

There is also a foundational motivation for this. Since 
the inception of quantum mechanics (QM), there has 
been an uneasy dichotomy between two points of view: 
is QM a fundamental description of nature or merely an 
algorithm to calculate probabilities for outcomes of ex- 
periments? The friction between these two viewpoints 
comes from the manifestly non-classical phenomena QM 
predicts, constrains or allows. This list of phenom- 
ena includes interference, the uncertainty principle, non- 
locality through entanglement (Bell inequalities) or oth- 
quantum telcportation and no-cloning theo- 
The relationships between these phenomena 
remain unclear Our interest lies in how the struc- 

ture of the allowed states of a quantum theory constrains 
its dynamics d, . We suggest that the results of this 
paper may be extrapolated to probabilistic theories more 
general than quantum mechanics. 

We utilize a symmetry of quantum states called envari- 
ance [Io| which emerges dynamically due to their entan- 
glement structure. This symmetry is a consequence of 
the tensor product structure (TPS) of quantum states^. 
In classical deterministic theories, Cartesian products are 



used to define assemblies of subsystems rather than ten- 
sor products and so a symmetry equivalent to envariance 
does not exist. 

We shall use an operational definition of the observ- 
ability of the change of a quantum state: the trace dis- 
tance (defined in Sec. V) between the transformed and 
untransformed state. An upper bound on this measure of 
the observability of the dynamics is derived. The bound 
is only dependent on the purity of the state. Another 
intriguing aspect of quantum states is the mathematical 
equivalence of states of subsystems of an entangled sys- 
tem to mixed states representing classical ensembles of 
quantum states. This connection allows the bound to 
apply to entangled systems as well. 

Envariance is defined in Sec. II. In Sec. Ill, we describe 
how the information contained in the subsystems of an 
entangled bipartite system can be less than the infor- 
mation contained in the whole system. It is also shown 
how the mixedness of a state constrains knowledge on 
all non-degenerate obscrvablcs of that state. These two 
qualities of quantum states are then used to motivate 
Sec. IV where the class of invariant operations on a com- 
pletely mixed state is considered. This symmetry is then 
used in Sec. V to derive an upper bound on the trace dis- 
tance between an untransformed state and a transformed 
one. Concluding remarks are presented in Sec. VI. 



II. ENVARIANCE 



* Electronic address: cacl.hassc@adclaidc.cdu.au 

1 It has been shown Q that TPSs are a generic feature of proba- 
bilistic theories with subsystems. Some form of envariance may 
then exist in such theories as well. 



Envariance [10|, LU| is a symmetry of entangled compos- 
ite systems. We define a composite system, in general, as 
a state that can be decomposed in terms of eigenstates 
of two or more mutually commuting observables where 
subsets of the total set of mutually commuting observ- 
ables completely describe subsystems. Therefore a parti- 
cle state with quantum numbers spin and position can be 
considered composite, with spin and position describing 
separate subsystems. 

Zurek's original use of envariance was to provide a 
proof of Born's rule under 'very mild' assumptions. We 
shall be using envariance in a different context and thus 
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assume Born's rule from the outset. 

Consider a composite system that can be decomposed 
into two subsystems a and /?, with the state \tp) £ H a <8 
Hp. Now suppose there exist unitary operators U a and 
Up, where 

U a :=U a ®lp, (1) 

with 1^ being the identity on Hp and U a : H a — > H a 
where H a is a subspace of H a , with an analogous defini- 
tion for Up. A state is said to be envariant under U a 
(or Up) if the following holds: 

U a Up\f) = \i>) or U a \i>) = Ul\f). (2) 

Note that from now on, we shall be considering finite 
dimensional Hilbert spaces only. 

Suppose we have a state of the form 

1 N 

IV>> = -^£ e ^K>l&>' (3) 

where {|oj)} and {|/3j)} form orthonormal bases for H a 
and Hp respectively, N = dim('H Q ,) = dim('H / 3) and <f>j 
are arbitrary phases. These states are envariant under 
all unitary transformations of H a (or Hp). When H a = 
H a , the state is maximally entangled (for subsystems a 
and P) and the group of envariant transformations is the 
group of all unitary transformations of H a alone (i.e., 
they can be decomposed as in Eq. 

Consider now the case of a state with Schmidt decom- 
position 

n 

i"> = E c ^K>ift>> ( 4 ) 

3=1 

with Cj G R + such that Cj ^ Cj for i ^ j, i.e., the coeffi- 
cients of the Schmidt decomposition have unequal norms. 
In this case, the group of envariant transformations in- 
cludes only relative (and overall) phase changes between 
the components |ai)|/3j), i.e., unitaries of the form 

n N 

U a :=(j2e i ^\a j )(a j \+ £ \a 3 )(a,\) ®lp, (5) 

i=l j=n+l 
n M 

Up:=(j2 e ~ iXk \fa)W+ E \Pk){Pk\)®Ic, (6) 

fe=l k=n+l 

where dim(H a ) = N,dim(Hp) = M and \j e (0,2tt) V 
j. These unitaries have the desirable property 

u a Up\n) = |0). (7) 

The most general case is where some of the coefficients 
Cj are equal and others are not. For the subspaces of H a 
spanned by the components whose coefficients are equal, 
we have envariance over the entire subspace. For the rest 
of the space, it is only relative phases of the components 
with unequal coefficients that can be envariantly trans- 
formed. 



III. ALLOWED STATES 

The emergence of envariance is a reflection of the prop- 
erty of entangled quantum states, where complete knowl- 
edge of the entire system (i.e., the state being pure) 
means incomplete knowledge of the subsystems. This 
can be understood in several ways: 

1. The reduced density matrices, tracing out a or /3 
(tr Q [p] or tr^p] for some pure p), have non-zero von Neu- 
mann entropy, leaving a mixed state partially equivalent 
to a classical lack of knowledge about the subsystem. 
However, the number of invariant degrees of freedom is 
only indirectly related to the amount of entanglement as 
mentioned in [llj . We shall consider this point in more 
detail later on. 

2. Consider a Bell state, 

W = ^(lt)lU>2 + U)l|t>2), (8) 

which is maximally entangled and as such has an SU(2) 
subgroup of envariant transformations which we can 
parametrize by the Pauli matrices, 

j l ^\v) = e il3 *\v), (9) 

where <r* is the i th Pauli matrix for the j th particle. Thus, 
rotating the spin of particle one is the same as rotat- 
ing the spin of particle two instead. This implies only 
the relative orientations of the rays within the subsys- 
tem Hilbert spaces are known. Zurck has cited a similar 
idea as his motivation for using envariance fl2l ] and calls 
it the 'relativity of quantum observables'. The situation 
can be said to have a kind of Machianity , analogous 
to the situation where the universe consists of point parti- 
cles and only relative distances between them are known, 
not global displacement or orientation. The state only 
contains information about the correlations between the 
particles. 

3. For maximally entangled subsystems, the probabil- 
ities for fine grained (non-degenerate) measurement out- 
comes of a subsystem whose reduced density matrix has 
its maximum von Neumann entropy become equal in any 
basis. This can be seen with the use of envariance, which 
is equivalent to a basis ambiguity of the subsystems. For 
instance, with the Bell state the probabilities for a par- 
ticular particle to be spin up or down in the z-direction 
are the same while the probabilities for the particle to 
be spin left or right in the x or y direction are also the 
same. This is in contrast to an uncntanglcd spin-I/2 par- 
ticle where there always exists a direction where the spin 
is definitely known. 

This last example is a special case of a phenomenon 
that does not apply to classical physics. The uncertainty 
principle is usually applied to pure states but the sit- 
uation changes for mixed states, such that the bounds 
on the uncertainties for incompatible observables become 
more strict. To show this, we utilize the concavity of the 
expression — xhix (with x G R + ) [3| such that for a den- 
sity matrix p and a fine grained basis = 5ij, 
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= the von Neumann entropy S(p) has the 

property: 

S(p) : = -tr[plnp] 

i 

<-J2(i\p\i)Hi\p\i)- (10) 

i 

Choosing {\i) \i} to be the eigenstates of a fine grained 
observable, then (i\p\i) is the probability to measure out- 
come "i" such that the Shannon entropy of said observ- 
able (call it O) is given by 

H (p) = -Y.(i\p\*)Mi\p\i)- (11) 

i 

Thus, 

S(p)<H (p). (12) 

For an alternative proof see [IH . 

This applies to all fine grained observables of the sys- 
tem described by p. In the case of the Bell state, the 
reduced density matrix obtained by tracing out one of 
the particles has maximum von Neumann entropy such 
that all non-trivial observables of the subsystem also have 
maximum Shannon entropy. Thus if a subsystem con- 
tains quantum correlations with another, the information 
we have about the subsystem is more constrained than in 
the case of classical physics where Shannon entropies of 
'incompatible' observables are allowed to be independent. 

IV. OBSERVABLE DYNAMICS OF 
COMPLETELY MIXED STATES 

Intuitively, when one lacks knowledge of a system, one 
expects our ability to distinguish the dynamics of the 
system to be lessened. We have seen that in the case of 
mixed quantum states, our knowledge of the system is 
less than allowed classically. 

We begin quantifying the distinguishability of dynam- 
ics of mixed states by extending envariance of completely 
mixed states to non-unitary operations. In this regard we 
choose to describe a quantum process in an operator-sum 
representation which maps density matrices to density 
matrices. A general physical operation on a can be de- 
scribed by a set [l6| of operation elements E^ a £ 'K a ®'K a 
where T-L a is the dual to H. The operation is then given 
by 

K 

£ a ( P ) = Y, E »°p E L (1 3 ) 

fe=l 

We shall be concerned with operations £ a whose ele- 
ments are trace preserving {^2 k E k E ka = I Q ) and also 
unital: 

K 

J2E ka El a =I a . (14) 
fc=i 



Let p a be a completely mixed state of a system a, purified 
by system /?, 

p a :=tr f3 W)m = ±I a . (15) 

All unital operations leave p a invariant (see Appendix 
A), i.e., 

£a{Pa) = Pa- (16) 

Some operations satisfying these conditions include: 

1. Unitary £ u {p) = UpW. 

2. Perfect premeasurements £p{p) = J2k=i PkpPk 
where Pk are projectors of a complete basis. 

3. Combinations of unitary and perfect premeasure- 
ment operations, e.g., £jjp{p) = £u £p(p) 

Interestingly, generalized measurements [l6j where the 
outcome is unknown do not necessarily satisfy these con- 
ditions, e.g.. for measurement operators of a two level 
system M 1 = |0)(0| and M 2 = |0)(1|, the left hand 
side of Eq. (|14[) with Ei a = Mj does not equal unity; 

m 1 m\ + m 2 m\ ^ I. 

V. UPPER BOUND FOR GENERAL MIXED 
STATES 

Let us now consider a general mixed state of a 

n 

p ol =Y J \c j \ 2 \a j ){a j \. (17) 
j'=i 

A purification of a by (3 is given by the Schmidt decom- 
position ([4]), where now the Cj and Cj may be equal for 
i =/= j. The group of envariant operations on |fi) is in 
general greatly reduced compared to the maximally en- 
tangled state. Thus the set of all operations that can be 
shown to leave p a invariant by the use of envariance is 
also reduced. This limits the previous proof of the unob- 
servability of the dynamics of a for cases where the state 
is not completely mixed. 

Our proposal is that even with a large reduction in the 
set of symmetries, the original set may apply in a partial 
sense. The motivation is that a large reduction in the 
symmetry can occur with only a very small reduction in 
the von Neumann entropy of a [ll| . Mixed states with 
less than maximum von Neumann entropy may still have 
some form of limitations on their observable dynamics 
for the full set of trace preserving operations satisfying 
Eq. (Tj"4")) . This turns out to be the case. 

To see this, we initially rewrite f2. Let us extend the 
su over j from one to N and define Cj = for n + 



2 However, the choice of the size of the extension may be chosen 
to be smaller depending on whether S leaves certain subspaces 
of Ha invariant. 
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1 < 3 < N. For M < N, we enlarge Up until the 
dimensionalities are equal. We then decompose into 
two parts, one that is maximally symmetric over W a and 
the rest of the state; 



Cj = — p= + dj 



(18) 



where dj :~ cj — l/y/N, such that 



N N 



Dchne 



1 N 



AT 



Q|fi 2 )=^e*K->|/3 ; 



(20) 
(21) 



where the constant Q = J^ZjdZ is chosen such that f2 



is normalized to 1. 

Our measure of the purity of a is given by Q. It is not 
equal to the usual measure of purity, which is tr[p^]. One 
can consider Q 2 as the x 2 value between the distribution 
of amplitudes {cj} and the constant distribution l/\/N. 
It follows from J2j c2 = 1 an d < cj < Cj that Q is 
bounded 



< Q < J2-2/VN. 



(22) 



The maximal value occurs for pure states of a, while 
Q = corresponds to completely mixed states (all Cj 
equal) . 

We now utilize a measure of the distinguishability of 
quantum states, the trace distance, defined as 



D(p,a) 



\tt\p-a\ 



(23) 



where p and a arc density matrices and \X\ := vXXX 
is the positive square root of X^X (defined by taking a 
spectral decomposition X^X = J^. ei\Xi)(xi\ and taking 

the positive square roots of the eigenvalues \f~XXX = 
It can be shown that UM 



D(p, a) — max tr[P(p — a) 



(24) 



where P is a projector and the maximization is taken 
over all possible projectors. This gives a clear physical 
interpretation of the trace distance. If an experimentalist 
wanted to distinguish whether they had the state p or a, 
the trace distance gives the maximum possible difference 
in probabilities for a projective measurement outcome for 
the two states. For instance, if for two states D = 1, it 
is in principle possible to do a projective measurement 



where the probability of getting a confirmatory result for 
one state is one while the other is zero and hence only 
one measurement is ever needed to distinguish the states. 

We are now in a position to derive an upper bound on 
D a := D(£ a (p a ),p a ). 

Let the state p n = |^)(^| be acted on by £ a as defined 
in Sec. IIVI In the case where Q = 0, |f2) = | Ox ) and 



£ a (p n |q=o) = £p (p n |q=o) , 



(25) 



wherc0 £p(p) = Yl<k=i E kPP E kf3 ( cf - Appendix A). As 
these two states are equal, they are indistinguishable. 
For the general case where Q may not be zero, a measure 
for the distinguishability of the two states can be given 
by the trace distance 



D 



Q/3 := D(s a (p a ),Sp(p n )y 



(26) 



In Appendix B, we show that D a p satisfies the following 
bound: 



D aP <2,/i-|i-g 2 



(27) 



This is related to D a in the following way. If f3 is an 
ancilla subsystem used to purify a or the experimentalist 
does not have access to subsystem /3, then we can ask 
about our ability to tell whether £ a has happened at all. 
This can be quantified by 

D(tvp[£ a (p n )},tr [£ (p n )}^ = D(£ a (tv p [p n }),tvp[p n }) 

= D a . (28) 

The partial trace over j3 is trace preserving, so D a is 
bounded by D a p: 



D a <D aP <2Jl- l-Q 2 + W 



(29) 



For values of Q < v2 — \/3 ~ 0.5, the right hand side 
of (|2T))) becomes less than one and hence bounds D a . 
For Q = 0, (|29p gives D a = which is the same result 
achieved in Sec. IIVI The upper bound on D a given by 
(|2"9l is our central result. 

The non-trivial nature of this bound can be seen by 
considering cases where D a is not bounded because Q is 
larger than s/2 — \/3: 

1. Consider two bases for a, {| a^)} and such 
that (a m \a m ) = for some m. Take a pure state 
c = \ot m )(otm\- The purity as given by Q is then, 



Q = J2-2/VN > V2-V3 



(30) 



3 For cases where "H p has to be enlarged and /3 is considered a real 
subsystem, Ep may not strictly be physical. 
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for N > 2. One can see that if a experiences a 
unitary transformation 



Otk\ 



(31) 



then the states have zero overlap: 

or D{U a aU a ,a) = 1. (32) 



\U a <rUi\ 



Thus the two states are in principle easily distin- 
guishable. 

2. Suppose the system a, still given by the pure state 
a, experiences a perfect premeasurement such that 



JV 



(33) 



where P; = \Ai)(Ai\ and {|^4i)} forms a com- 
plete orthonormal basis for H a such that \a m ) = 

SfcLn(VV^V)|-A»)- It is convenient to use the def- 
inition of fidelity [l6| for density matrices p and 

T 



F(p, t) := tr [VpVVpV^] , (34) 

to obtain 

(P^^PaP)) 2 = (« m |^P l( rP|a m ) 



1 

TV' 



(35) 



In this case, the fidelity bounds the trace distance 
l-(P(a,^PaP)) 2 



l-±<D(*,Y,Pi<7Pi)- 



(36) 



Thus, the observability of the process 
D{cr, PivPi) tends to 1 as N tends to oo. 

Finally, we note that the bound may be extended to 
mixed states of a composite a and (3 system. Let 



E 



(37) 



where y r TO = 1 and p m is a pure density matrix of 
the composite system Vm. Define Q m as the Q measure 
of the purity of the tvp(p m ) state and define £ a and £p 
in the usual way. Then, using the convexity of the trace 
distance, 



D\£ a (p) 1 £p(p)j < r m D [£ a (p m ), £p (p m ) 

m 

< 2 £r 



(38) 



Thus, the distinguishability of the dynamics of a in a 
mixed composite state is bounded by the average of the 
bounds of the pure states p m . 



VI. REMARKS 

We have shown that given a trace preserving unital 
operation, the trace distance between the transformed 
state and its original is bounded by given that the 
purity is enough (Q < 0.5). For maximally mixed states 
where Q = 0, the bound implies the operation must be 
unobservable. 

The bound is motivated in Scc.QTDand Sec.lTVlon 
the intuition that lack of knowledge of a state leads to 
lack of an ability to distinguish the dynamics. We note 
that trace preserving, unital operations cannot decrease 
the von Neumann entropy [17j . This leads us to ask 
whether the class of trace preserving, unital operations 
is the largest such class where (|2l?|) or a stronger bound 
holds that depends only on the purity of the input state. 
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Appendix A: Symmetries of Completely Mixed 
States 



The proof of the invariance of p a under unital oper- 
ations is trivial. Here we provide an alternative proof 
which give the tools needed for Sec. |V] The first step is 
to extend the symmetry of the second version of ([2]) for 
pure state ((3|), to Pa]^) = A t l s|'0) where p a := p a (8)1/3 is 
a general linear operation on pure states. The operations 
p a could for instance be a projector onto a subspace of 
T~i a . Also suppose \tp) is maximally entangled with re- 
spect to subsystems a and /?. Since p a acts identically 
on subsystem /3, it follows that 



/v 



3=1 



(39) 



Define {oti\p a \otj) := pij such that 



N N 
3=1 i=l 



N N 



V i=l j=l 

1 N 



(40) 
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where :— J2f=i VijlPj) ■ For non- unitary p a , 

{\/3j}} need not be orthonormalQ 
Define 



N 



Hp :=Ia®^e #fc |/3fc)(/?fc| := la <8> p-p, (41) 



fe=l 



(42) 



This symmetry is not equivalent to envariance as pp may 
not be invertible and hence in general there does not exist 
a [ip 1 such that p, a p,^ \ip) = \ijj). With Eq. (pf2"j) . we can 
consider symmetry properties of completely mixed states. 
Let the completely mixed state to be a state of system 
a; 



Pa :=txpU){ip\] = ^l ol . 



(43) 



Suppose we have a quantum operation £ a (p) = p' that is 
given in an operator-sum representation; 



K 



(44) 



k=l 



where E ka are linear maps E ka : T-L a — > 1-L a . The effect 
of £ a upon p a is then 



K 



£M = J2E ka trp[\ij){^\}El a 

fc=l 

K 

= trp[Y d E ka m^\El 



(45) 



fc=i 



Since £?fc a is of the form there exists an E k p such 

that EfcalV') = E k p\ip) and similarly (-01 -^L = ^\ E lp- 
Then 



if 

^a(Pa)=trp[5^£? W |V>(V|^ 
fc=l 
K 

= tTp[Y,ElpE k p\*l>)^\ 

k=l 



^Y: k= i E lpE k p=lp,then 



£ a {Pa) = Pa- 



(46) 



(47) 



The completely mixed state is then invariant under £ a 
if £ a satisfies the condition J^fcLi E\oE k p = Ip, which is 
equivalent to the unital condition 



A" 



^2 E ka El a - I Q 



(48) 



k=l 



Many \f}j) may even be the null state. 



Appendix B: Proof of upper bound of trace distance 

To prove Eq. (|2"T1) . we shall need a few identities. Using 
the definition of d; and taking to be normalized, we 
find; 



N N . 



N 



t=l 
N 



■di 



N 



JV 



dj 



i=i 

We also compute the overlap 
using Eq. (|49| . Defining 



(49) 



AT 



(50) 



Pll 


:= |fii)(fii 


P12 


:= |ni)(fia 


P21 


:= |fi 2 )(Q 1 


P22 


:= |n 2 )<fi2 



(51) 



we have 

p = |n><n| 

= Pll + Q(Pl2 + P2l) + Q 2 P22- 

The ps have some useful relationships, namely, 



(52) 



2 2 

Pll = PW1P22 = P22, 

P12P21 = p\\,P\2pYl = - 

P21P12 = P22,P2lP21 = - 



Q 

-fPl2, 
"IP- 



(53) 



The trace distance has particular properties which shall 
also be used. Thus, 

D a p: = D(£ a (p),£p{p)^ 

< D(£ a (p),£ a (p n )^ +D(£ a (p n ),£p(pj) 

= D {e a (p) , £ a ( Pll )) + D (£p (pn ) , £p (p)) 

< D(p,pn)+D( pi i,p) 



= 2£>(p u ,p). 
Because p and pn are pure, 



D 



(pn,p) = \Jl - (F{Pn,p)) 2 , 



(54) 



(55) 
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where F is the fidelity; 

7-1/ n , \( 1/2 l/2\l/2 



Thus, 



,V2 



(pil[Pll + Q(Pl2 + P2l) + Q 2 P22]P11 

(pu - \Q z pw - ^Q z ph 



1/2 



D aP <2Jl- 1-Q 2 



kQ 4 \ 



(57) 



l n 2 n i/n2„ , ln4„ \V 2 



l-Q 2 



I04 



(56) 



[1] I. L. Chuang and M. A. Nielsen, J. Mod. Optics 44, 2455 
(1997). 



[2] M. A. Nielsen, eprint |arXiv:quant-ph/9606012| (1996). 
[3] A. Gilchrist, N. K. Langford and M. A. Nielsen, Phys. 

Rev. A 71, 062310 (2005). 
[4] C. H. Bennett, D. P. DiVincenzo, C. A. Fuchs, T. Mor, 

E. Rains, P. W. Shor, J. A. Smolin and W. K. Wootters, 

Phys. Rev. A 59, 1070 (1999). 
[5] W. K. Wootters and W. H. Zurek, Nature, 299, 802 

(1982). 

[6] J. Oppenheim and S. Wehner, Science 330, 1072 (2010). 
[7] S. Popescu and D. Rohrlich, Found. Phys. 24, 379 (1994). 
[8] J. Barrett, Phys. Rev. A 75, 032304 (2007). 



[9] 

[10] 

[11] 
[12] 



[13] 
[14] 
[15] 
[16] 



[17] 



D. Gross, M. Miiller, R. Colbeck and O. C. O. Dahlsten, 

Phys. Rev. Lett. 104, 080402 (2010). 

W. H. Zurek, Phys. Rev. A 71, 052105 (2005). 

M. G. A. Paris, J. Quant. Inf. 3, 655 (2005). 

J. A. Wheeler and W. H. Zurek, Quantum Theory and 

Measurement, p. 772, (Princeton, USA: University Press, 

1983, Princeton Series In Physics). 

J. B. Barbour, Class. Quant. Grav. 11, 2853 (1994). 

A. Wehrl, Rev. Mod. Phys. 50, 221 (1978). 

B. W. Schumacher, Phys. Rev. A 54, 2614 (1996). 

M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information, (Cambridge University Press, 
Cambridge, 2004). 

D. Perez-Garcia, M. M. Wolf, D. Petz and M. B. Ruskai, 
J. Math. Phys. 47, 083506 (2006). 



